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Abstract. The complete magnetic and multiferroic phase diagram of Mni_ x Co x W04 
' O ' single crystals is investigated by means of magnetic, heat capacity, and polarization 

experiments. We show that the ferroelectric polarization ~P in the multiferroic state 
q | abruptly changes its direction twice upon increasing Co content, x. At x c i=0.075, 

P rotates from the b— axis into the a — c plane and at x C 2=0.15 it flips back to the 
. b— axis. The origin of the multiple polarization flops is identified as an effect of the 

^ ' Co anisotropy on the orientation and shape of the spin helix leading to thermody- 

namic instabilities caused by the decrease of the magnitude of the polarization in 
, the corresponding phases. A qualitative description of the ferroelectric polarization 

\^ ' is derived by taking into account the intrachain (c— axis) as well as the interchain 

■ (a— axis) exchange pathways connecting the magnetic ions. In a narrow Co concen- 

tration range (0.1<x<0.15), an intermediate phase, sandwiched between the collinear 
high-temperature and the helical low-temperature phases, is discovered. The new phase 
. . . exhibits a collinear and commensurate spin modulation similar to the low-temperature 

^ 1 magnetic structure of MnWO.4. 
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1. Introduction 

The coexistence of different fundamental orders in solids and their mutual interaction 
have attracted attention of physicists and chemists for more than a century. 
Magnetoelectricity and the coupling between magnetism and ferroelectricity have 
inspired experimentalists and theorists because of many novel physica phenomena that 
arise from the complex cross correlations between magnetic and electric orders as 
well as the potential of these multiferroic materials for future applications in form of 
magnetoelectric sensors or new memory elements [1, 2, 3]. 

With the discovery of ferroelectricity induced as a secondary order parameter 
by an inversion symmetry breaking magnetic order in TbMn03 [4], the interest in 
multiferroic magnetoelectric compounds has been revived and the research in the field 
has made significant progress since then. The transverse helical spin structure is 
one of the candidates to induce ferroelectricity whenever the spin lattice coupling is 
strong enough to cause the cooperative ionic displacements resulting in a macroscopic 
electrical polarization. Some typical examples of spin-helical multiferroics are TbMn0 3 
[5], Ni,3V208 [6], and MnWC>4 [7, 8]. The important role of helical magnetic orders 
in inducing ferroelectricity in multiferroics has been theoretically studied within the 
Ginzburg-Landau model [5, 9] or microscopic theories [10, 11]. 

The origin of helical magnetic orders is frequently found in magnetic frustration, 
either due to geometric constraints (for example antiferromagnetic interactions on a 
triangular lattice) or because of competing exchange interactions. In frustrated systems, 
there are several magnetically ordered phases that are very close in energy and they 
compete for the ground state. Small perturbations can therefore control the magnetic 
phases and stabilize one or the other of the competing orders. This results in an extreme 
sensitivity of most multiferroic materials with respect to magnetic or electric fields 
[8, 12, 13, 14, 15], external pressure [16, 17, 18, 19], and ionic substitutions [20, 21]. 
While the magnetic field couples to the ordered moments affecting their orientation 
and changing the magnetic structure, the application of external pressure compresses 
the lattice and tunes the microscopic exchange interactions and anisotropy parameters. 
Whereas the control of microscopic interactions via the application of pressure is indirect 
and limited by the lattice compressibility, the substitution of the magnetic ions by other 
ions with different spins, exchange couplings, and anisotropy constants allows for a more 
targeted control of the magnetic interactions. 

The ideal multiferroic compound for investigating the effects of substitutions is 
MnWO"4. Also known as the mineral Hiibnerite, MnWO"4 crystallizes in a monoclinic 
structure (space group P 2/c). At low temperature, several magnetic phase transitions 
and multiferroic properties have been reported [7, 8]. The spins of the only magnetic 
ion, Mn 2+ , form different frustrated orders upon decreasing temperature because of 
competing magnetic exchange interactions. A sinusoidal magnetic order with an 
incommensurate modulation defined by the vector C^ 3 = (—0.214,1/2,0.457) sets in 
at Ttv=13.5 K [22]. The Mn-spins of this AF3 phase are collinear and they are confined 
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to the a — c plane at an angle of 34° with the a— axis. At Tf=12.6 K, the spins form 
a non-collinear order by tilting partially toward the b— axis. The magnetic order in 
this AF2 phase is described by a spin-helix with the same propagation vector, £^ 2 = ^3- 
Since the helical order breaks the spatial inversion symmetry, the AF2 phase is also 
ferroelectric with a polarization vector pointing along the 6-axis, in perfect agreement 
with the spin-current model [10] and the Ginzburg-Landau theory [9]. Only at lower 
temperature, at T L =7.5 K, the magnetic order locks into a commensurate structure 
described by the vector = (-0.25,0.5,0.5). The order in this AF1 phase shows the 
highly frustrated tt44 s P m sequence which preserves the spatial inversion symmetry. 
Therefore, the AF1 phase is paraelectric. 

The complex phase sequence is a direct consequence of a high degree of frustration 
in the magnetic subsystem. The frustration originates from competing exchange 
interactions [22] and results in several magnetic states having almost the same energy 
and competing for the ground state. Recent neutron scattering experiments have 
revealed the long-range character of the magnetic interactions involving up to 11 different 
exchange coupling constants and the spin anisotropy to explain the low-energy magnetic 
excitations [23]. Replacing the magnetic Mn ion by other transition metal ions with 
different spin, exchange coupling, and anisotropy parameters will tune the various 
magnetic and multiferroic phases and, possibly, reveal new physical phenomena in 
frustrated systems. Mn WO 4 is a good candidate since it forms stable compounds when 
Mn is completely replaced by Fe [24], Co, Ni [25], or Zn [26]. A uniform solid solution 
of MnW04 with (Fe,Co,Ni,Zn)W04 is expected to form over the whole concentration 
range because all ternary compounds are isostructural with relatively small variations 
of the lattice constants. 

The substitution of Fe 2+ (spin 2) for Mn 2+ (spin 5/2) has shown an interesting 
magnetic phase diagram [27] with coexisting phases at some critical Fe concentrations 
[28]. The multiferroic properties of Mni_ x Fe x W0 4 have been investigated earlier [29] 
and it was shown that the substitution of more than 4 % Fe suppresses the multiferroic 
AF2 phase. On the contrary, the replacement of nonmagnetic Zn 2+ in Mni_ x Zn r W04 
did stabilize the AF2 phase as the ground state and the commensurate (tt44-) AF1 
phase was completely suppressed [30, 31]. A similar effect was also observed in the 
Co substituted compound, Mni_ :r Co x W04, based on magnetic and neutron scattering 
experiments of polycrystalline samples [32]. The neutron scattering data further 
suggested the rotation of the FE polarization from the b— axis toward the a— axis at 
higher Co content. A sizable a— axis component of the polarization was indeed recently 
found in a single crystal of Mno.9Coo.1WO4 [33, 34]. At slightly higher doping (15 % Co), 
however, the only component of the FE polarization was aligned with the b— axis [35]. 
These seemingly conflicting results indicate an extreme complexity of the multiferroic 
phase diagram of Mni__ r Co ;c W04 and warrant further exploration of single crystals of 
various Co concentrations. 

In order to study the multiferroic and magnetic states and the polarization rotations 
between b— and a— axes of Mn 1 _, E Co :E W04 in more detail, we have investigated the full 
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doping range of Mn 1 _ a .Co x W04 between x=0 and x=0.3 through magnetic, dielectric 
(polarization), and heat capacity measurements. We show that, at zero magnetic field 
and upon increasing Co concentration, the ferroelectric polarization rotates twice, from 
b into the a — c plane at x~ 0.075 and back to b for x> 0.135, limiting the multiferroic 
phase with the polarization perpendicular to b to a narrow region in the phase diagram. 

2. Experimental 

A large number of single crystals of Mni_- r Co ;E W04 with fourteen different compositions 
x between and 0.3 have been grown in a floating zone optical furnace. Powder X-ray 
diffraction was used to check the phase purity of the polycrystalline feed rod before 
the crystal growth. No impurity phases could be detected within the resolution of the 
spectra. The chemical composition and the Co content of the single crystals were verified 
by energy-dispersive X-ray (EDX) measurements testing up to 15 different spots of a 
single crystal. For the 2 % Co sample the actual Co content was 2.6 %, slightly higher 
than the nominal substitution level. For all other samples the Co content was close to the 
nominal composition within the error limits of the EDX data. In the following sections, 
we will use the nominal composition to distinguish between different substitution levels. 

The crystals cleave easily and expose the a — c plane which was used as the contact 
area for b— axis polarization measurements. For measurements along other directions 
(a— and c— axes), the crystals were cut accordingly. Single crystal Laue diffraction was 
employed to determine and confirm the crystalline orientations. 

For dielectric and ferroelectric polarization measurements, small plate-like samples 
were cleaved or cut from the bigger single crystal. Silver paint was used for electrical 
contacts. The contact area was about 10 to 15 mm 2 and the sample thickness 
was typically less than 0.5 mm. The sample was mounted on a home-made probe 
inserted into the Physical Property Measurement System (Quantum Design) for precise 
control of temperature and the speed of temperature change. The spontaneous 
polarization was measured through the pyroelectric current method using the K6517A 
electrometer (Keithley). The experimental protocol was chosen to diminish any artificial 
contributions to the current as follows: The sample was cooled from 30 K (above Ttv) 
in an electric field of about 3 kV/cm to low temperature while measuring the current. 
After releasing the applied voltage and shortening the contacts for several minutes at the 
lowest temperature, the pyroelectric current was measured upon increasing temperature. 
The rate of temperature change was 1 K/min in all experiments. The polarization is 
than calculated by integrating the current from high temperature (in the paraelectric 
state) to low temperatures. The cooling data for the polarization have been found 
consistent with the (zero electric field) warming data, but they reveal large thermal 
hysteresis effects in some part of the phase diagram. 

The bias electric field during cooling (3 kV/cm) has to be high enough to ensure 
the alignment of the FE domains. This was verified by measuring the FE polarization 
after cooling in a reduced electric field of 2 kV/cm resulting in the same magnitude of 



5 



the FE polarization. Some experiments had been conducted in cooling in a negative 
(sign-reversed) electric field and the obtained polarization was of equal magnitude but 
of opposite sign, as expected in a FE state. 

Magnetization measurements have been conducted in a commercial magnetometer 
(Magnetic Property Measurement System, Quantum Design). The low-field data were 
typically taken in a magnetic field of 200 Oe. The heat capacity was measured using 
the relaxation method implemented in the heat capacity option of the PPMS. 

3. Results and Discussion 

3.1. Magnetic susceptibility of Mn\- x Co x WO± 

The small amount of Co doping of only 0.02 suppresses the low-temperature paraelectric 
AF1 phase. This is revealed in magnetization data shown in Fig. la. The transitions 
into the sinusoidal (Tjv) and the helical magnetic phases (T<^) are resolved in anomalies 
of the b— axis magnetic susceptibility, as indicated in Fig. la. The phase transition from 
the helical AF2 to the commensurate AF1 phase (Tl) is only detected in MnW0 4 as a 
sharp increase of Xfe(T) at T L . In Mn 98 Co .02WO 4 , however, this transition is already 
missing. The quick suppression of the AF1 phase extends the ferroelectric AF2 phase 
to the lowest temperature. This result is consistent with similar conclusions derived 
from the study of polycrystalline Mni_ x Co x W0 4 [32]. An analogous doping effect on 
the stability of the AF1 and AF2 phases has been observed recently in single crystals 
of Mn 1 _ a; Zn ;E W0 4 [31]. For 0.02<x< 0.05, Xb(V changes very little (different curves in 
Fig. la are vertically offset) but T^ and Tc are slightly decreasing. 

With increasing x between 0.075 and 0.15, however, Xb(T) changes significantly 
below Tc- The b— axis magnetic susceptibilities in this range of Co substitutions are 
shown in Fig. lb. The relatively sharp decrease of Xb below T c disappears above 
x=0.075 and it changes into a plateau-like feature that dominates the low-temperature 
behavior of the susceptibility. The low-T plateau exhibits additional anomalies in form of 
small step-like changes of Xfc(T) for x>0.12 suggesting additional transitions or changes 
of the magnetic orders and a more complex magnetic phase diagram in this doping range 
(details are discussed below). 

With further increasing x>0.17, however, the T-dependence of Xb again changes 
qualitatively, as shown in Fig. lc. While the Neel transition is clearly defined by a sharp 
slope change of x&(T), the susceptibility continues to increase to lower T until it suddenly 
decreases sharply below a second transition temperature, T c . The T-dependence of Xb 
below T c is similar to the susceptibility data shown in Fig. la for the low-doping range. 
The two anomalies labeled as Tjv and T c separate two magnetic phases. With the 
results from powder neutron scattering [32] and our earlier results for x=0.15 [35], the 
phase between T^ and Tc can be identified as the AF4 phase with the commensurate 
magnetic modulation corresponding to (^ 4 = (0.5,0,0), similar to the magnetic phase 
observed in Fe-substituted MnW0 4 above 12.5 % iron content [27]. 
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Figure 1. 6— axis magnetic susceptibility of Mni_ x Co x W04. The arrows indicate 
the critical temperatures Tjv (AF3, AF4), T c (AF2), and T L (AF1). Different curves 
are vertically offset for better clarity by 0.35, 0.3, and 0.4 units in (a), (b), and (c), 
respectively. 



3.2. Heat capacity of Mn\- x Go x WO± 

The various phase transitions between different magnetic phases are reflected in distinct 
peaks or anomalies of the heat capacity, C P (T). For example, the heat capacity of the 
undoped MnWC>4 exhibits pronounced peaks at all three magnetic transitions from the 
paramagnetic phase to AF3, from AF3 to AF2, and from AF2 to AF1 [7]. Our single 
crystal data for x=0 are included in Fig. 2a, top curve. For 0.02<x< 0.1, the heat 
capacity data, shown in Fig. 2a, reveal only two peaks between 11 and 14 K, representing 
the transitions from the paramagnetic to the sinusoidal and subsequently to the helical 
phases, respectively, similar to MnWC>4 (note that the AF1 phase is already suppressed 
in this doping range). 

Upon increasing x, the onset of magnetic order at TV is accompanied by a step-like 
increase of C p and a large peak of C p develops near 11 K together with a small third 
peak at lower temperature between x=0.12 and x=0.15 (Fig. 2b). For x>0.17, the low- 
temperature peak of C P (T) disappears and only two anomalies remain, indicating that 
two magnetic phases dominate the phase diagram in this doping range (see Fig. 2b). 
The two transition temperatures exhibit a large x-dependence with the higher transition 
(Tjv) quickly increasing and the second transition (Tc) shifting to lower temperatures. 

The x- dependencies of the low-temperature properties of Xb and C p suggest that 
there are at least three doping ranges (denoted by I, II, and III) in the phase diagram 
of Mn 1 _ a ,Co a; W04 between x=0 and x=0.3 that are distinguished by their characteristic 
magnetic properties and thermodynamic anomalies, as shown in Figs. 1 and 2, 
respectively. In ranges I (0.02<x<0.05) and III (0.17<x<0.3) only two transitions 
are observed, separating two phases with different magnetic structures. In range II, 
the phase sequence upon decreasing temperature appears to be more complex involving 
more than two magnetic orders. Therefore the question arises which of the various 
phases is ferroelectric (or multiferroic). This can only be decided by measurements of 
the ferroelectric polarization. 
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Figure 2. Molar heat capacity of Mni-^Co^WC^. Different curves are vertically 
offset for better clarity by -0.5 units with reference to the x=0 and x=0.12 data in (a) 
and (b), respectively. The labels in (a) for x=0 indicate the peaks associated with the 
magnetic transitions at T N (AF3), T c (AF2), and T L (AF1). 



3. 3. Ferroelectricity in Mri\_ x Co x WO4 

The ferroelectric polarization i*(T) in the AF2 phase of MnWC>4 is directed along the 
monoclinic b— axis [7]. This is in accordance with the symmetry and the prediction of 
the spin current model, with the b— axis lying in the plane of the spin helix and assuming 
that the major contribution arises from the average nearest neighbor spin-spin exchange 
oriented along the c— axis. For x=0.1, however, it was shown that the polarization 
had rotated by 90° into the a — c plane [33]. A smaller c— axis component found in 
Mno.9Coo.1WO4 was attributed to the magnetic interaction along the a— axis between 
spins of different chains [33, 34]. In order to account for all possible orientations of i*(T) 
for 0<x<0.3, we have measured for each composition the polarization components along 
a, b, and c. 

For Co substitution levels up to 0.05, ^(T) was found to be aligned with the b— axis 
(Fig. 3a). The undoped MnW04 shows the sharp drop of to zero at the lock-in 

transition into the AF1 phase, in agreement with earlier reports [7, 8]. With x increasing 
above 0.02, however, the paraelectric AF1 phase is quenched and the ferroelectric AF2 
phase extends to the ground state. The low-temperature (measured at 5 K) polarization 
decreases with x and reaches about 50 % of the value of MnW04 at x=0.05. The 
reference value for MnW04 at 5 K, Pj, ~60 /zC/m 2 , was derived from data measured 
at a magnetic field of 3 T which also suppresses the paramagnetic AF1 phase [7]. The 
ferroelectric critical temperature Tc decreased slightly from 12.8 K (x=0) to 12.3 K 
(x=0.05). At the upper limit of this Co concentration interval, for x=0.042 and 0.05, a 
small component of the polarization along the a— axis was detected. The size of P a was 
more than one order of magnitude smaller than Pf,. This small P a could be an effect of 
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Figure 3. Ferroelectric polarization of Mn 1 _ a .Co x W04. (a) Low doping range 
(x<0.05, P b ), (b, c) intermediate doping (0.075<x<0.135, P a and P c ), and (d) high 
doping range (x>0.15, P&). For x=0.12, warming and cooling data are presented 
to show the thermal hysteresis. The dashed line for x=0.075 shows the maximum 
polarization achieved in magnetic field. 



a minute misalignment of the crystal or it could indicate the onset of the polarization 
flop from b into the a — c plane. Our measurements of the ferroelectric polarization 
of Mn .95Coo.o5W0 4 in magnetic fields [36] are consistent with the onset of a smooth 
rotation of the polarization (and the spin helix) at doping levels as low as 0.042. 

At slightly higher Co doping, from x=0.075 to x=0.135, the ferroelectric 
polarization does not show any significant b— axis component within the resolution of 
the measurement. Small values of detected in some measurements can be due to 
a minute misalignment (typically less than 3°) from the perfect orientation. However, 
a large polarization was found along the a— axis accompanied by a smaller component 
along c, as shown in Figs. 3b and 3c, respectively. The magnitude of P c is too large to 
be explained by a misalignment and it has to be understood as an intrinsic effect. Since 
the nearest neighbor intrachain magnetic couplings are oriented along the c— axis, they 
cannot generate a polarization along c and the interchain magnetic interactions along a 
have to be considered. This is discussed in more detail in Section 3.5. 

The magnitudes of P a and P c at 5 K decrease quickly with increasing x and nearly 
vanish at x=0.135. For x=0.075, both P a and P c pass through a maximum at 9 K and 
decrease upon further cooling. This unusual behavior can be attributed to a temperature 
dependent reorientation of the helical plane from the optimal position in the a — c plane 
resulting in a decrease of the polarization below 9 K. A magnetic field oriented along 
the b— axis was shown to stabilize the a— axis component of (as for example in 

Mn .95Coo.o5W04 [36]). Indeed, P a of Mn .925Coo.o75W04 measured in a field H b =5 
Tesla rises continuously and reaches the maximum value of more than 120 /iC/m 2 at 
5 K, as shown by the dashed line in Fig. 3b. It should be noted that this value is 
the highest one measured in any doped or undoped MnW04 compound and it exceeds 
the maximum b— axis polarization of MnW04 by a factor of two. It is also interesting 
that this maximum of ~P appears right at the critical x at which the rotation of the 
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Figure 4. Low-temperature (5 K) values of the ferroelectric polarization of 
Mni-zCozWC^ showing the separation into three ferroelectric phases I, II, and III, 
with different orientations of P , For x=0, the value of P was extrapolated to low 
temperature based on data in magnetic fields suppressing the AF1 phase. For x=0.135 
the maximum value obtained at 5 Tesla is shown as a star. 




polarization from the b— axis into the a — c plane happens. 

The data for x=0.1 are in quantitative agreement with recent reports [33, 34]. 
Upon further increasing x, P(T) decreases quickly and the critical temperature of the 
ferroelectric phase shows a sharp decrease from 11.3 K (x=0.1) to 7.6 K (x=0.135). A 
significant thermal hysteresis of more than 1 K develops above x=0.1, as demonstrated 
by the warming and cooling data for x=0.12, included in Fig. 3b. The values of Tq 
of the ferroelectric phase match perfectly with the temperatures of the small peak of 
the heat capacity (Fig. 2b). Near the upper limit of region II of the phase diagram 
(x=0.135), the polarization at 5 K dropped already to below 2 fiC/m 2 . The fast decrease 
of ~P with increasing x and the fact that no sizable component of P was observed along 
b indicate a major change of the spin helix in a narrow doping range. 

Increasing the Co content x above 0.135 results in another flop of the ferroelectric 
polarization from the a — c plane back to the b— axis with a sudden increase of its 
magnitude. Data for x>0.15 are shown in Fig. 3d. The results for x=0.15, previously 
reported by us [35], are included as the dashed line. The recovery of the 6-axis 
polarization up to 30 /iC/m 2 at 5 K suggests another sudden rotation of the spin helix 
in Mni_ x Co x W04 at a critical doping of 0.15. No sizable a— or c— axis components of 
P could be detected in this higher doping range. 

The x-dependency of the ferroelectric polarization is summarized in Fig. 4 showing 
the values of 7^ at 5 K vs. Co doping x. The two instabilities near x cl =0.075 and 
x c2 =0.135 are clearly represented and it is obvious that the rotation of P happens 
when its magnitude drops below a critical value. Above x=0.3, ~~P decreases quickly to 
zero and the helical magnetic phase is quenched completely. The collinear AF4 phase 
with the characteristic spin modulation defined by = (0.5, 0, 0) becomes the ground 
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Figure 5. Magnetic and multifcrroic phase diagram of Mni_ x Co x W04. Phase 
boundaries are derived from magnetic, heat capacity, and polarization measurements. 
The preliminary phase assignments are based on recent neutron scattering reports [38] . 
The shaded area indicates the intermediate paraelectric AF1 phase, as discussed in the 
text. 



state. This is consistent with the AF4 magnetic structure that is realized for x=l, i.e. 
in C0WO4 [37]. 

3. 4- Magnetic and multiferroic phase diagram of Mni- x Co x W04 

The polarization measurements discussed above define the ferroelectric (or multiferroic) 
phases of Mni_ a; Co a ;W04. However, other magnetic phases, e.g. the paraelectric AF1, 
AF3, or AF4 phases, have to be defined through their characteristic properties and the 
magnetic as well as thermodynamic anomalies at the phase boundaries (see sections 
3.1 and 3.2 above). Combining all experimental data, we can completely resolve the 
magnetic and multiferroic phase diagram as shown in Fig. 5. The phase assignment of 
different ordered structures used in Fig. 5 is based on previous [32] and most recent [38] 
neutron scattering results. 

Some features of the phase diagram, such as the quick suppression of the AF1 phase 
and the overall onset of magnetic order at T^r, are in agreement with the results from 
powder neutron scattering experiments [32]. The Neel temperature decreases slightly 
with x for x<0.1, the first magnetic structure in this doping range is the sinusoidal AF3 
phase. However, T^ increases significantly with x>0.1 when the commensurate AF4 
phase replaces the incommensurate AF3 phase. The transition from the sinusoidal AF3 
to the helical AF2 phase with the polarization along the b— axis extends to x~0.075. 
Between x c i=0.075 and x C 2=0.135, however, the ferroelectric polarization is constrained 
to the a — c plane exclusively (AF5 phase). 

It should also be noted that for 0.12<x<0.15 three magnetic phase transitions 
are observed in magnetic data (Fig. 1) and heat capacity (Fig. 2). A third phase 
(shaded area in Fig. 5), not reported so far, is sandwiched between the AF4 and AF5 



11 



phases. According to the polarization measurements, this phase is paraelectric. The 
transition into this phase is accompanied by the large peak of the heat capacity, as 
shown in Fig. 2b. Recent detailed neutron scattering experiments [38] have revealed 
the magnetic structure of this phase, identifying it as a commensurate collinear magnetic 
structure with the wave vector of the low-temperature AF1 phase of MnWOzi. This is 
also consistent with our previous neutron data for x=0.15 [35]. The details, however, 
have yet to be explored. The presence of this intermediate phase also explains the sudden 
decrease of the ferroelectric transition temperature in this range and its recovery above 
x=0.15. 

With further increasing x (>0.15), only two phase transitions can be observed. At 
Ttv, a commensurate magnetic order (AF4 phase) develops first, as reflected in the kink 
of Xb (Fig. lc) and the high-temperature peak of C p (Fig. 2b). The second transition 
at Tc recovers the b— axis polarization, the magnetic structure therefore has to be 
noncollinear allowing for a polarization along b exclusively, similar to the AF2 phase 
at lower doping levels. as well as the critical temperature T c decrease quickly with 
x to P 6 (2 K)=0.6 /iC/m 2 and T c =3.6 K, respectively, at x=0.3. The strong reduction 
of the polarization may be due to a decrease of the magnetic moments contributing to 
the spin helix, a possible change of the ellipticity, or to the increasing influence of the 
AF4 magnetic structure. A possible coexistence of both phases was suggested based on 
powder neutron scattering experiments [32]. 

3.5. Discussion and possible origin of the multiple polarization flops in Mni_ x Co x WO4 

The multiferroic state in MnWC>4 is the classical example of a polar structure induced by 
an inversion symmetry breaking helical magnetic order [7, 9, 15]. The strong coupling of 
the lattice to the spin helix, as revealed e.g. through significant anomalies of the thermal 
expansivity [39], results in the polar distortion of the lattice and the ferroelectric state. 
Both, the microscopic theory as well as symmetry arguments relate the magnitude and 
the direction of the electrical polarization to the vector product of the term (rjtj x rrlj) 
with the position vector i^- connecting both magnetic moments at sites i and j along 
the propagation direction of the spin helix. Taking into account the specific structure 
of MnWC>4 with two spins per unit cell, the magnetic moment in the helical state can 
be expressed as: [22] 

rflil^iu) = rft\\cos(!f~Fl i + $ Q ) + r^ ± sin(~^^i + $ Q ) (1) 

itj Q is the position vector to the Mn site a (=1, 2) in the unit cell i, if is the 
propagation vector of the spin helix defined as g« = 27r/Aj (Aj is the wavelenght of the 
spin modulation along the i th direction), and $ 2 = ( £i+<2 , zc/2-|-7r. The two perpendicular 
vectors rfty and rft± are chosen to define the helical magnetization plane with my and 
mi being equal to the long and short half axis of the magnetic ellipse, respectively. 

To derive a general expression for the ferroelectric polarization we have to calculate 
the vector product of the two magnetization vectors at the pair of atomic sites coupled 
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by the magnetic exchange interaction. The result for nearest neighbor spins within a 
chain along the c— axis is: 

r?f (itj 2 ) x rr£(i^ji) = m\\m±sin(q z c/2)lt (2) 

ft is the normal vector of the helical plane (length 1), as indicated in Fig. 6. 

As noted above, the coupling between spins exclusively within the chains 
(c— axis) cannot account for the observed c— axis component of the polarization in 
the intermediate doping range. Magnetic exchange between different chains has to 
contribute to the ferroelectric polarization. This is also suggested by the results of 
recent inelastic neutron scattering data [23] which lead to an estimate of the Heisenberg 
exchange constants for different pairs of spins along as well as in between the chains. 
The magnetic exchange between spins along the b— direction was generally very weak; 
however, sizable exchange constants had been determined for neighboring spins along 
the c— as well as a— axes. Therefore, it appears conceivable to include contributions to 
the polarization that arise from the spin-spin interactions along the a— axis: 

rrt(^ i+ i ja ) x rft(l$i a ) = m\\m±sin(q x a)jf (3) 

Here itj+i, a and if iQ refer to two magnetic ions in neighboring chains displaced 
along the a— axis. 

The contributions to the electrical polarization can now be expressed as the vector 
product of the expressions (2) and (3) with the corresponding average position vectors 
connecting the spins in each case. Using spherical coordinates, we express the normal 
vector it = sin(Q)cos(ip)e x + sin(Q)sin((p)e y + cos(Q)e z . is the angle of ~n with 
the z-axis and <p denotes the angle of the x-y plane projection of it with the x-axis. 
it i are the unit vectors along the i-th cartesian coordinate (i=x,y,z). For simplicity, we 
will consider ~~£ ' x parallel to it and lt z parallel to it (the angle between it and it is 
nearly 90°), with ~~£ ' y parallel to b . The polarization in terms of 6 and (p due to the 
intrachain exchange is: 

j£( ) _ (j^m\\m±sin(q z c/2)[—sinQsinip~^' x + sinQcostplt y ] (4) 

A similar expression is obtained for the interchain exchange along the a-axis: 

v (2) 

p =C^m\\m±sin(q x a)[—cosOlty + sinOsin(p~~(t z \ (5) 

and are constants not depending on 6 or (p with supposedly |C^ 2 ^| < |C^^|, 
because of the larger distance between the chains and the weaker magnetic exchange. 
The signs of q z and q x determine the helicity of the spiral spin modulation. 

In MnWC>4 the helical plane is perpendicular to the a — c plane forming an angle 
of about 34° with the a— axis, corresponding to the spin easy axis in the collinear AF1 
and AF3 phases [22]. Accordingly, the orientation of the helix is defined by 0=34° and 
(p=180°. This results in a b— axis polarization with contributions from equation (4) as 
well as equation (5). Whether the two contributions to Pt are additive or subtractive 
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(a) . .(b). 

c Pocnxc\\b C Pxnxc\\a 




0<x<0.5 0.075<x<0.135 

Figure 6. Schematic diagram of the possible orientations (defined by the normal 
vector #) of the spin helix in Mn1_a.Coa.WO4 for x<0.05 (a) and 0.075<x<0.135 (b). 
The shaded area in (a) symbolizes the range where the helical plane comes close to the 
a— axis and the associated b— axis polarization is too small to stabilize the AF2 phase. 



depends on the helicities of the spin spiral and the signs of the prefactors in equations 
(4) and (5), i.e. the microscopic details of the magnetic interactions. 

The dependence of the magnitude of the ferroelectric polarization, on the angle 
6 was discussed by Sagayama et al. [15] taking into account only the contribution 
due to nearest neighbors along the c— axis, i.e. equation (4). However, because of the 
second contribution from equation (5), the angular dependence of P&(6, ip) will be more 
complex. 

In the following we will qualitatively discuss the complex x-dependent ferroelectric 
polarization as a consequence of the possible changes of the spin helix induced by 
the cobalt substitution. Using formulas (4) and (5), our focus will be on the angular 
dependence, i.e. the expressions in the square brackets, although the possible changes 
of the prefactors with x may cause additional changes of the magnitude of ~P. The 
reorientation of the helical plane was recently observed for Co concentrations between 
0.05 and 0.1 [33, 34]. The driving mechanism, however, has yet to be revealed. 

It appears conceivable that the change of the average spin anisotropy upon doping 
with cobalt plays a major role. The easy axis of the transition metal (T) spin in TWO4 
strongly depends on the T-ion [40] . For C0WO4 it was determined that the angle of the 
easy axis with the a— axis is 6 = —46° [37], a change of 80° with respect to MnW04. 
It can therefore be expected that the Co substitution tends to rotate the easy axis and 
the helical plane in the multiferroic phase of Mni_ a; Co 2 ;W04 into the direction of the 
spin easy axis of C0WO4, as schematically indicated by the curved arrow in Fig. 6. 
This rotation was indeed observed in neutron scattering experiments for small x<0.05 
[32, 33, 38]. The decrease of from 34° (x=0) to 15° (x=0.05) results in a reduction 
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of the ferroelectric polarization without changing its orientation, according to the data 
shown in Fig. 3a and Fig. 4. 

The decrease of P&, however, is energetically unfavorable since the multiferroic state 
gains energy from the coupling of the polarization to the magnetic order parameter. 
The energy gain originates from a third order term in the free energy expansion that 
is proportional to the polarization. This term, which has to be invariant with respect 
to the spatial inversion and time reversal operations, has the mathematical form of a 
Lifshitz invariant: 7^ • [aI(V V)aI] [9]. The gain of free energy in the 

multiferroic state scales accordingly with the magnitude of the polarization. If the spin 
helix rotates closer to the a— axis (6 <15°), the polarization as well as the energy gain 
in the multiferroic state will further decrease as long as the angle p remains unchanged. 
For energy reasons, we can therefore expect an instability at a concentration x c i when G 
decreases below a critical value, C , and the multiferroic AF2 state cannot gain enough 
energy to be stable with respect to other competing states. The shaded area in Fig. 6a 
shows schematically the region of instability for the helical plane in this phase. 

The multiferroic state can gain more energy through a flop of the helical plane 
allowing for larger values of the ferroelectric polarization. The flop can be visualized as 
a rotation of the helical plane toward the a — c plane (for (p=90° and 6=90° the helical 
plane coincides with the a — c plane). According to equations (4) and (5), this results 
in a rotation of P into the a — c plane, where the two components P a and P c arise from 
equations (4) and (5), respectively. It is interesting to note that the angular dependence 
of both components is exactly the same implying a similar temperature or magnetic 
field dependence with different magnitudes of the polarization values (according to the 
different prefactors in equations (4) and (5)). 

The a — c polarized state is realized in the AF5 phase between x cl =0.075 and 
x c2 =0.135 (Figs. 3b and 3c). The first signature of the polarization flop in the phase 
diagram of Mni^Co^WC^ appears at x=0.075 (Fig. 3b). P a assumes the largest values 
above 8 K, exceeding even the values for x=0.1 (Fig. 3b and previous reports [33, 34]). 
With further cooling, however, P a passes through a maximum and decreases again. This 
unusual behavior may indicate a temperature dependent tilt of the helical plane below 
8 K away from the optimal orientation. The c— axis component exactly mimics the 
temperature dependence of P a at a smaller scale, as is expected from equations (4) and 
(5). Due to this low-temperature tilt, the maximum polarization cannot be obtained 
for x=0.075 in zero magnetic field. A magnetic field applied along the b— axis, however, 
reveals the maximum value of P a >120 fiC/m 2 at 4.5 K, which appears to be the upper 
limit for the MnW0 4 system. The high-field data are included in Fig. 3b (dashed line) 
and Fig. 4 (star). 

P a and P c decrease quickly between x cl and x c2 and no sizable polarization 
component along the b— axis was observed. The data for x=0.1 are in good agreement 
with previous reports [33, 34]. The observed decrease of the polarization may indicate 
another rotation of the spin helix, so that 6 <90° and/or ip <90° and sinQsimp < 1, 
as sketched schematically in Fig. 6b . However, the decrease of P a and P c to almost 
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zero (P a =2/iC/m 2 for x=0.135 and T=5 K) without any sizable b— axis polarization 
is difficult to understand. If ©, ip systematically decrease with x, both equations 
(4) and (5) are expected to also give rise to a significant component P&, which was 
not experimentally observed. The question of whether the two contributions to P5 
accidentally may cancel each other has yet to be explored. Another possibility for 
the decrease of P(x) in this section of the phase diagram could be an increase of the 
ellipticity (my m±) of the spin helix with x, reducing the prefactor in equations (4) 
and (5). If m±=0, the magnetic order is collinear and no polarization is allowed. The 
true nature of the magnetic order in the AF5 phase can only be revealed by detailed 
neutron scattering experiments. 

Nevertheless, the decrease of P a and P c again diminish the stabilizing energy of 
the multiferroic state arising from the magnetoelectric coupling term in the free energy 
expansion. Therefore, another instability causes the polarization to flop back toward the 
6— axis (for x>0.15). A sizable value of the ferroelectric polarization along the b— axis is 
recovered and the multiferroic state is stabilized with regard to the nearly paraelectric 
state close to x=0.135. The further decrease of P b is driven by the competition with 
the commensurate and paraelectric AF4 magnetic structure that extends to lower 
temperatures and coexists with the helical phase below T c . The phase boundary is right 
at x=0.15, as previous neutron scattering data of Mn a85 Co .i5WO4 [35] have revealed 
the coexistence of a sizable amount of AF5 phase with the arising AF2 phase (with 
b— axis polarization). From the above discussion it is now obvious that the AF5 phase 
at x=0.15 (which was shown to be the major phase at low temperatures) does not 
generate any significant polarization, but it reduces the AF2 phase to a minor phase. 
This explains the small value of P;, (shown in Fig. 3c) as compared to the values for 
x>0.17. 

4. Summary 

We have investigated the complete multiferroic and magnetic phase diagram of 
Mni_ a; Co a ;W04 for 0<x<0.3 through magnetization, heat capacity, and polarization 
measurements. With increasing x, the ferroelectric polarization abruptly changes 
direction at the critical Co concentrations of x cl =0.075 (P flips from the 6— axis into 
the a — c plane) and x c2 =0.15 (P flips back from the a — c plane to the 6— axis). Above 
x=0.3, the multiferroic state ceases to exist and the only magnetically ordered state 
is the commensurate and collinear AF4 phase, the same phase as in C0WO4. The 
change of magnitude and orientation of the ferroelectric polarization as function of x is 
explained by the rotation of the spin helix driven by the anisotropy of the magnetic ions 
(Co 2+ vs. Mn 2+ ). The driving mechanism of the multiple polarization flops is suggested 
to be of thermodynamic origin and it is related to the reduction of the magnetoelectric 
energy gain, whenever the rotation of the spin helix causes a significant decrease of 
the polarization. A model relating the magnitude and direction of the polarization to 
the orientation of the spin helix is derived and provides a qualitative description of the 
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multiferroic properties in the phase diagram. In particular, it is shown that the nearest 
neighbor intrachain magnetic coupling alone cannot reproduce the experimental results 
and the longer range interchain interactions have to be taken into account. Between 
x=0.1 and x=0.15, a new magnetic and paraelectric phase is discovered, sandwiched 
between the AF4 and AF5 phases. The existence of this phase is the origin of the 
observed dip in the T^(x) function of the multiferroic phase. 

Note added: After completion of this work we have been informed about a related 
work by Urcelay-Olabarria et al. [41] on Mno.sCoo.2WO4 (x=0.2). The reported results 
for the magnetization and ferroelectric polarization are in good agreement with the data 
reported here. 
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